Interaction of colloids with a nematic-isotropic interface 
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The Landau-de Gennes free energy is used to calculate the interaction between long cylindrical 
colloids and the nematic-isotropic (NI) interface. This interaction has two contributions: one is 
specific of liquid crystals and results from the deformation of the director field close to the particles 
or to the interface, while the other is generic and results from wetting and surface tension effects. 

Deep in the nematic phase the director field of long cylindrical colloids, with strong homeotropic 
anchoring, exhibits two half-integer defect lines. As the colloid moves towards the interface, the 
director configuration changes through a series of discontinuous transitions, where one or two of the 
defects are annihilated. In addition, the NI interface bends towards the colloid in order to minimize 
the elastic free energy in the nematic. In the isotropic phase, the colloid is surrounded by a thin 
nematic layer that reduces the surface free energy under favorable wetting conditions. 

The interaction has a well-defined minimum near the interface. In this region the director and 
interfacial structures are complex and cannot be described analytically. Using the numerical results 
for the Landau-de Gennes free energy in the harmonic region, we obtained simple scaling laws for 
the (linear) force on the colloid. 

PACS numbers: 61.30.Cz, 61.30.Jf, 61.20.Ja, 07.05.Tp 



I. INTRODUCTION 

Colloidal dispersions are suspensions of solid or liquid 
particles in a host fluid. The size of the particles can vary 
from a few nm up to several /im. Colloidal dispersions are 
often found in long-lived metastable states, providing the 
basis for a range of industrial applications, e. g. paints, 
drugs, foods, coatings, etc [|. 

Colloidal dispersions in nematic liquid crystals form 
a special class of colloids. The difference from ordinary 
colloids arises from the long-range orientational order of 
the liquid crystal molecules, described by the so-called 
nematic director. Topological defects of the director 
field 0, additional long-range forces between colloidal 
particles Jp, 0, 0] , and, as a result, supermolecular struc- 
tures [Ulij are phenomena specific of colloidal nematics. 

Colloidal nematics are usually prepared in the isotropic 
phase. To prevent flocculation due to attractive 
van der Waals forces, colloidal particles are treated to 
induce electrostatic or steric repulsive interactions. Af- 
ter cooling to a temperature below the NI transition, the 
colloids often segregate forming non- uniform clusters. In 
general, topological defects stabilize the dispersion [j| al- 
though, in some cases, they may help flocculation by pro- 
ducing additional attractive interactions [3. 

The spatial distribution of colloidal particles is very 
sensitive to the cooling rate and cooling conditions 0, |j . 
Recent results [T(| revealed that the drag on colloids by 
a NI interface plays an important role in the final spatial 
pattern. In the experiment, large nematic and isotropic 
domains were separated by moving NI interfaces. De- 
pending on the cooling rate (velocity of the interface) 
different structures were observed: a) cellular structures, 



also reported earlier 0, 0, ^| , with particle- free nematic 
domains separated by particle-rich regions; b) striped 
structures, where the particle-rich regions are arranged 
in a set of stripes; c) root-like structures. A major con- 
clusion of the study was that the motion of the particles 
may be controlled by their interaction with the NI inter- 
face. 

Here we reexamine this conclusion by analyzing in 
detail the interaction between colloidal particles and a 
NI interface. As a first step we investigate the (static) 
Landau-de Gennes order parameter distributions around 
colloidal particles in the vicinity of NI interfaces. By 
calculating the corresponding free energies we obtain ex- 
plicit results for the forces. 

For simplicity, we consider the interaction of a long 
cylindrical colloid with a planar NI interface. We assume 
strong homeotropic anchoring of the director at the parti- 
cle's surface everywhere and wetting boundary conditions 
(see section III C|l . In the nematic phase, these anchoring 
conditions give rise to a pair of 1/2 strength defect lines, 
that merge with the isotropic phase as the particle ap- 
proaches the NI interface, distorting the planar interfacial 
region on a scale of the order of the colloidal radius. As a 
result, the interface is strongly bent around the colloidal 
particle on both sides of the interface. The effective in- 
teraction is shown to be rather complex and to exhibit 
a well defined minimum close to the interface. Near the 
minimum, a scaling analysis of the numerical results for 
the free energy yields the scaling form of the force on the 
particle. 

The paper is organized as follows. The geometry, the 
model, and the minimization technique are described in 
Sec. [H] The results for the director configurations and 
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FIG. 1: xy cross-section, of size L x L, of the system studied. 

free energy are presented in Sec. IIIII A discussion of the 
major findings, together with concluding remarks, appear 
in Sec. El 

II. LANDAU-DE GENNES THEORY 
A. Geometry 



where the direction of maximal orientational order is 
given by the director n, Q is the scalar order parame- 
ter, and the unit vectors l,n,m form a local orthonormal 
triad. 

Spatial non-uniformities that do not coincide with n 
break the cylindrical symmetry of the average angular 
environment, and thus one must allow for biaxiality, i.e. 
B ^ 0. 

The translational symmetry of the interface along the 
z axis and the assumed homeotropic boundary conditions 
on the colloidal surface, imply that the director is con- 
fined to the xy plane. Then the vector m may be chosen 
along the z axis and the vector I is in the xy plane 

n = (cos 6, sin 6,0), (2) 
I = (sint9, — cost9, 0), 
m = (0,0,1). 

In this case, the tensor order parameter has three in- 
dependent components only 

f Qu Qu o \ 
Q = Q12 Q22 (3) 
\ -Qn-Qaa / 

In the following we will use this representation of the 
tensor order parameter. 



In this study we consider a geometry similar to that 
used in the experiments 0|- Briefly, one end of the 
sample is placed in a hot and the other in a cold oven. 
The temperatures are chosen such that the NI interface is 
in the gap between the ovens and a constant temperature 
gradient is maintained throughout the sample. 

A geometry mimicking the experimental setup is shown 
in Fig. ^ A uniform temperature gradient is imposed 
along the x axis with the 'hot' wall at x = L/2 and the 
'cold' wall at x = —L/2, where L is the size of the system. 
The director at the 'cold' wall is fixed either parallel or 
perpendicular to the wall. The order parameter at the 
'hot' wall is set to zero. Note, that the position of the in- 
terface, at x — 0, is pinned by the temperature gradient. 

A colloidal particle, which we take to be a long cylin- 
der of radius R, with the symmetry axis parallel to the 
z axis, is immersed in the nematic. Periodic boundary 
conditions are used in the y direction. 

B. Tensor order parameter 

We use the Landau-de Gennes tensor order-parameter 
Q a /3 formalism since it is free from divergent terms due 
to the defect cores and it takes into account the biaxiality 
that may occur in non- uniform nematics [l2l IT3I Il4| . 

Owing to the traceless, symmetric character of the ten- 
sor order parameter it can be represented as 0] 

Q a p{r) = ^Q{3n a np - 8 a p) + -B(ljp - rriamp), (1) 



C Free energy 

The system is described by the Landau-de Gennes free 
energy [l^l 

HQ} =f(fb + fe)dV + J fsdS (4) 

where fb is the bulk free energy density, f e is the elas- 
tic free energy density and f s is the surface free en- 
ergy. Within a mesoscopic approach the minimum of 
the Landau-de Gennes functional T gives the equilibrium 
value of the tensor order parameter. 

Symmetry arg uments yield for the local bulk free en- 
ergy density [Mll3 

f b = aTrQ 2 - WrQ 3 + c [TrQ 2 ] 2 , (5) 

where a is assumed to depend linearly on the tempera- 
ture, while the positive constants b, c are taken tempera- 
ture independent. 

It is convenient to scale out the variables by defining 

Qij = 6c/bQij, (6) 

f b = 24 2 c 3 /6 4 / 6 - 

It will be understood that such scaling has been carried 
out, and we shall omit the overbars in the text below. 

We also introduce a dimcnsionlcss temperature r by 
defining 

a = rb 2 /2Ac. (7) 
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For a uniform uniaxial nematic (Qn — Q, Q22 = 
Q33 = —1/2Q) the free energy 10 takes the form 

h = tQ 2 - 2Q 3 + Q 4 . (8) 

The nematic state is stable when r < 1 with a degree of 
orientational order given by 




We modeled the temperature gradient in the x direction 
by assuming that r depends on the x coordinate 

T = r c (l + a^j , (10) 

where r c = 1 is the NI transition temperature. Equa- 
tion (|10|l implies that the NI transition occurs at x = 0, 
i.e., the interface is in the middle of the cell and is parallel 
to the y axis. 

The elastic free energy density can be written as 01 

2 dxk dxk 2 dxj dxk ' 

where the constants L\ and L2 are related to Frank- 
Oseen elastic constants by K\\ = = 9Q 2 (Li + 

L^/2)/2 and K22 = QQ 2 L\/2 and Qb is the bulk nematic 
order parameter. The sign of L2 defines the preferred 
orientation of the director at the NI interface. L2 > 
(L2 < 0) favors planar (perpendicular) anchoring [l8| . 

We assumed strong homeotropic anchoring of the di- 
rector at the colloidal surface. This is valid if the anchor- 
ing parameter WR/K » 1, where W is the anchor- 
ing energy of the surface, and holds for large colloidal 
particles and/or anchoring strengths. We also assumed 
that the nematic at the particle surface is uniaxial with 
a scalar order parameter Q s = 1, independent of x. Un- 
der these conditions, a planar interface is wetted by the 
nematic phase as the temperature is lowered through the 
NI transition. 



D. Minimization procedure 

The equilibrium distribution of the tensor order pa- 
rameter Qij is obtained by minimizing the free energy 
functional J2J numerically using finite elements with 
adaptive meshes. During the minimization the square in- 
tegration region L x L was triangulated using the BL2D 
subroutine 19] . The functions Qij are set at all vertices 
of the mesh and are linearly interpolated within each tri- 
angle. The free energy is then minimized using the conju- 
gate gradients method [2(j under the constraints imposed 
by the boundary conditions. 

A new adapted mesh is generated iteratively from the 
previous minimization. The new local triangle sizes are 
calculated from the variations of the free energy, in or- 
der to guarantee a constant numerical weight for each 




1 > 1 1 1 1 1 1 1 1 1 

-0.2 -0.1 0.1 0.2 0.3 

x 



FIG. 2: Order parameter Q and biaxiality B of a free NI 
interface. Circles: order parameter, results of the numerical 
minimization; solid line: fit to Eq. 1121 . Squares: biaxiality, 
numerical results. Note the different scales and offset for the 
order parameter and biaxiality profiles. 



minimization variable. The final meshes, with a mini- 
mal length of ~ 10~ 2 , had ~ 10 5 minimization variables. 
Lengths are in units of the diameter of the colloidal par- 
ticle, 2R=1. 

We used £ 2 = L^c/b 2 = 0.001 that corresponds to a 
nematic correlation length £ w 0.03, which is larger than 
the smallest mesh length. The system size was either 
L = 12 or L = 20 and the temperature gradient was set 
at a = 0.025. 

In order to obtain both stable and metastable con- 
figurations, we used different types of initial conditions 
including flat and curved NI interfaces. 

III. RESULTS 
A. Free interface 

First, we analyze the tensor order parameter of a free 
NI interface, without colloidal particles and temperature 
gradients, for a system with L2/L1 = 2. For this ra- 
tio of the elastic constants the interface favors planar 
anchoringof the director (molecules parallel to the in- 
terface) [lg ■ The order parameter is fixed at the cell 
boundaries: Q(x = —L/2) = 1, Q(x = L/2) = and 
the director is aligned along the y axis, parallel to the 
interface. 

The order parameter and biaxiality profiles are shown 
in Fig. |21 As expected, the order parameter changes from 
the bulk value of the nematic phase, Qb = 1, to the bulk 
value of the isotropic phase, Qiso = 0. The width of 
the interface is of the order of the nematic correlation 
length. The nematic and isotropic phases are uniaxial, 
with zero biaxiality. The interfacial region, however, is 
slightly biaxial. 

Note, that the analytic solution of the correspond- 
ing variational problem does not exist, even in the one- 



dimensional case, when the tensor order parameter de- 
pends only on the x coordinate. However, one may use 
the de Gennes Ansatz for the order parameter profile, 
which neglects biaxiality [Isj . 

Q = i(l-tanh^, (12) 

where £ is the nematic correlation length when the direc- 
tor is parallel to the interface 

C 2 = e (e + . (is) 

A fit of the results of the numerical minimization to 
Eq. J32J) yields ( 2 ~ 0.00807 while Eq. JT3J gives C 2 = 
0.008, i.e., the agreement is very good. The small differ- 
ence is attributed to the biaxiality of the interface ^3 > 
not taken into account by the Ansatz 1|12[) . Note that 
we decreased the minimal mesh length down to 10~ 3 in 
order to increase the accuracy of the results. 

We also found that the surface tension is a function of 
the polar angle of the director with the interface normal. 
As a result, the interface has an easy axis with corre- 
sponding anchoring energy. We found planar interfacial 
anchoring (easy axis parallel to the interface) for > 0. 
L-2 < results in homeotropic interfacial anchoring (easy 
axis normal to the interface) in agreement with previous 
results [Hill. 

We studied the interaction of a colloid with the NI in- 
terface under various anchoring conditions of the direc- 
tor at the interface: i) director parallel to the interface, 
L2 > 0; ii) perpendicular to the interface, L2 < 0; iii) 
director tilted with respect to the interface, due to a mis- 
match between interfacial anchoring and the alignment 
at the cell boundaries. 



B. Planar interfacial anchoring 




FIG. 3: Order parameter and director maps for colloids at a 
distance, d, from the NI interface: (a) d = —2.5; (b) d = —2; 
(c) d = -1; (d) d = 1.5; (e) d = 1.625; (f) d = 2. Red 
(dark) corresponds to the isotropic phase with Q = 0, and 
white to the nematic phase with Q«l. System size L = 12, 
anisotropy of the elastic constants L2/L1 = 2. 



The anisotropy of the elastic constants is LijL\ = 2, 
favoring director alignment parallel to the NI interface. 
The director at the 'cold' wall is fixed parallel the y axis, 
which is also parallel to the NI interface. 

Typical order parameter and director maps are shown 
in Fig. The strong anchoring of the director at the 
colloidal surface yields two half-integer defect lines in the 
nematic, far from the interface, to ensure that the topo- 
logical charge of the system is zero (see Fig. |3^). The 
director distortion vanishes very rapidly in the nematic 
phase, and the cores of the defects extend over a few 
nematic correlation lengths, in line with previous stud- 
ies 0, El E3 • At these distances there is no interaction 
between the colloid and the interface: the director at the 
interface is uniform, the interface is flat and it is pinned 
at x — 0, where the nematic and isotropic phases coexist 

( T = 1 )- 

On reducing the colloidal distance from the interface, 
the defect closest to the interface merges (discontinu- 
ously) with the isotropic phase (Fig.|2l3,c). The interface 



bulges towards the colloid to accommodate the isotropic 
phase, where the defect core disappeared. Note that the 
anchoring at the NI interface is planar and follows the 
interfacial curvature, except in the region where the de- 
fect used to be. Here the director tilts and the interfacial 
anchoring becomes homeotropic, see Fig. [3Jd. 

As the colloid moves further into the isotropic phase, it 
is wrapped by the NI interface that forms a nematic 'cav- 
ity' around the colloid (Fig. 01). The second defect is still 
present on the nematic side. At a certain point, this con- 
figuration becomes metastable, and eventually annihila- 
tion of the second defect occurs. This is accompanied by 
a symmetry breaking transition: the cavity is no longer 
symmetric under y — > —y reflexion (Fig. 03 illustrates one 
of the two possible configurations). These configurations 
with broken symmetry are always metastable. 

Finally, deep in the isotropic phase, the colloid is 
wrapped by a thin layer of nematic phase due to the 
(wetting) boundary conditions at the colloidal surface 
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FIG. 4: Free energy T as a function of the distance of a 
colloidal particle from the NI interface. Both metastable and 
stable solutions are shown. Circles: particle moving from 
the isotropic to the nematic phase; squares: particle moving 
from the nematic to the isotropic phase. System size L = 
12, anisotropy of the elastic constants L2/L1 = 2. Letters 
correspond to the director configurations shown in Fig. |3] 



(Q s = 1), Fig.Ef. 

To investigate the nature of the structural transitions 
between different director field configurations, we plot 
the free energy as a function of the particle position in 
Fig. H 

Deep in the isotropic phase, the configuration of Fig-Ef 
is the only stable one. There is a small force acting on 
the particle in the direction opposite to the temperature 
gradient. This force is due to an increase in the free 
energy of the nematic layer, as the particle moves into a 
region at higher temperature. 

Moving from the isotropic into the nematic phase, this 
configuration becomes first metastable, and then unsta- 
ble. At a certain point, the configuration with a nematic 
'cavity' in the isotropic phase, Fig. I2JI, becomes stable. 

Finally in the nematic phase, the configuration with 
one defect becomes metastable, and then unstable. The 
second defect appears, as shown in Fig. Oil- At this point 
the interaction of the particle with interface almost van- 
ishes. The small negative slope of the free energy, giving 
rise to a force along the temperature gradient, is due to 
the temperature dependence of the elastic energy, that 
increases as the particle moves towards the cold wall. 

The jumps in the free energy reveal discontinuous 
(or first-order) instabilities between the physical director 
configurations. The equilibrium position, i.e., the posi- 
tion where the force on the colloid vanishes, is near the 
interface, at x = 0. 



C. Homeotropic interfacial anchoring 



Liquid crystals with negative 
anisotropy, K 2 2 > Ku,K; 
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FIG. 5: Order parameter and director maps for colloids at 
a distance, d, from the NI interface: (a,b,c) d — —1.5; (e) 
d = —0.5. System size L — 20, anisotropy of the elastic 
constants L2/L1 = —0.5. 



[l6| . In these systems the NI interface favors homeotropic 
anchoring, i.e., alignment of the molecules perpendicular 
to the interface [l]| . 

To study this case, we fixed the director at the cold 
wall along the x axis and used for the elastic constants 
L2/L1 = —0.5. For this set of parameters the director is 
parallel to the x axis, or perpendicular to the interface. 

The order parameter and director maps are shown in 
Fig. [3 Since the director is perpendicular to the in- 
terface, the defect lines are in the plane parallel to the 
interface, Fig. [S^. In this geometry both defects merge 
with the isotropic phase very quickly, when the parti- 
cle is still on the nematic side of the interface, Fig. [SJ3. 
Subsequently, the particle is surrounded by a thin layer 
of nematic phase until the fully wrapped state in the 
isotropic phase obtains, Fig. |3f. 

Deep in the nematic phase, we have also found a config- 
uration with only one defect, Fig. 03. However, analysis 
of the free energy, Fig. shows that this configuration 
is always metastable. 

D. Director tilted at the NI interface 

Another situation, often found in experiments, corre- 
sponds to tilted NI interfaces. However, the (simple) 
Landau-de Gennes theory fails to describe this interfacial 
anchoring and higher order gradients or more sophisti- 
cated density-functional theories are required in order to 
describe tilt at NI interfaces [23|. In finite cells, on the 
other hand, tilted interfaces may result from mismatch- 
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FIG. 6: Free energy T as a function of the distance of a 
colloidal particle from the NI interface. Squares: initial con- 
dition with a 'fiat' interface; circles: initial condition with a 
'curved' interface; triangles: 'curved' interface with the initial 
director tilted at 7r/4 (metastable configuration). System size 
L — 20, anisotropy of the elastic constants L2/L1 = —0.5. 
Letters correspond to the director configurations shown in 
Fig. 



ing boundary conditions at the cell boundaries and the NI 
interface proper. Indeed, in a planar cell, one may have 
the director aligned parallel to the bounding plates, say, 
along the x axis and a NI interface in the zy plane pro- 
viding planar anchorin g along the z axis. This is the so- 
called hybrid, or 7r-cell |2J, |25j where the director bends 
in the bulk to match the boundary conditions; as a re- 
sult there is a bulk splay-bend deformation even in the 
absence of colloidal particles. 

For this geometry we used L2/L1 = 2 and fixed the 
director at the 'cold' wall parallel to the x axis or per- 
pendicular to the NI interface. The system size is L = 20. 

The order parameter and director maps, shown in 
Fig. [7| are qualitatively similar to those described in sec- 
tion ^^Bj where the director is parallel to the cold wall 
and the interface. The tilt of the defects is due to the 
tilted nematic director, Fig. EJi. Otherwise, the same 
configurations are present: with two defects, deep in the 
nematic phase (Fig. UJi); with one defect, close to the 
interface (Fig. Et>,d); without defects, also close to the 
interface (Fig. Ek); without defects, deep in the isotropic 
phase (similar to the configuration shown in Fig. |3f ) • 

However, a careful analysis of the free energy, which is 
shown in Fig. [SI reveals that the structural transitions are 
different. The configuration with one defect (Fig. 01) be- 
comes metastable in the isotropic side of the interface; the 
defect- free configuration (Fig. (7^), that was metastable 
previously, is now stable. A physical explanation is sim- 
ple: because of the tilt, the remaining defect is closer 
to the interface and it is more easily annihilated. The 
bending of the interface is also weaker. 
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FIG. 7: Order parameter and director maps for colloids at a 
distance, d, from the NI interface: (a) d = —2.5; (b) d = —2; 
(d,e). Red (dark) corresponds to the isotropic phase with 
Q — and white to the nematic phase with Q ~ 1. System 
size L = 20, anisotropy of the elastic constants L2/L1 = 2. 
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FIG. 8: Free energy T as a function of the distance of a 
colloidal particle from the NI interface. Squares: initial con- 
dition with a 'flat' interface; circles: initial condition with 
a 'curved' interface. System size L = 20, anisotropy of the 
elastic constant L2/L1 — 2. The director is tilted at the NI 
interface. Letters correspond to the director configurations 
shown in Fig. [7| 



IV. DISCUSSION AND CONCLUSIONS 

We end with a brief discussion of the implications of 
our work on the drag of colloidal particles by the NI 
interface. 

The free energy, Figs. 01 and reveals that close 



7 




-1 -0.5 0.5 
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FIG. 9: Force on the particle as a function of the particle 
position. Inset: strength of the force as a function of the 
colloidal size and fit to k oc R 2 . 

to the minimum T is a quadratic function of the posi- 
tion of the colloid and thus, in this region, the force is 
proportional to the separation, d. 

We have calculated the force F = —dF/dd, as a func- 
tion of d, for several colloidal sizes, in the geometry with 
planar interfacial anchoring, as in the section ITlI Bl The 
resulting curves, together with the fit to the following 
dependence 

F = Fo + k li (W) 

are shown in Fig. |5J The inset of Fig. [5] shows the de- 
pendence of the force strength k on the particle size. We 
find that, to a good approximation, k oc R 2 . 

It is interesting to compare this result to the scaling 
arguments used in Ref. [lfj . If we neglect phenomena re- 
lated to the creation and annihilation of defects, and treat 
the interface as sharp, the only parameters in the problem 
are the radius and length of the colloid, R and L respec- 
tively, the elastic constant K, and the surface tension a. 
The appropriate combinations of these parameters with 
the dimensions of force are KL/R and aLR, where we 
took into account that for a particle of length L >> R 
the force is proportional to L. Therefore, the elastic and 



the surface tension contributions to the total force can be 
written as F c = KL/Rf e (d/R), F s = Laf s (d/R), where 
/ e and f s are dimensionless functions of the penetration 
depth, d/R. It is clear that neither of these forces gives 
the observed oc R 2 dependence on the particle size: the 
elastic force is inversely proportional to R, while the con- 
tribution due to the surface tension does not depend on 
the particle size at all. One reason for this discrepancy 
(apart from the presence of defects) is the complex struc- 
ture of the interface, that bends forming a cavity to wrap 
the colloid. 

Another interesting conclusion is the asymmetry of 
the free energy profile as a function of the colloid po- 
sition: when the particle moves from the nematic to the 
isotropic phase, the metastable phase with a nematic 
'cavity' (Fig. 01) becomes unstable much faster than the 
corresponding metastable state with one or no defects 
(Figs. Eb), when the particle moves from the isotropic 
to the nematic phase. Therefore, the interface drags the 
particle more efficiently when the colloid moves from the 
nematic to the isotropic side, owing to the presence of 
a higher energy barrier. This effect has been observed 
experimentally |Toj . 

We also note that wrapping a colloidal particle by the 
NI interface resembles the process of wrapping a colloidal 
particle by a membrane [2g, [23 ■ A detailed comparison 
of these processes requires, however, further analysis. 

To summarize, we studied the interaction of a cylin- 
drical colloidal particle with a NI interface. We found 
that, when the particle is close to the interface, the force 
on the particle is proportional to the penetration depth 
d/R, and the amplitude of the force scales as R 2 . At 
larger penetrations, discontinuous configurational transi- 
tions occur, related to the creation/annihilation of topo- 
logical defects. 
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